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Abstract 

We have studied the rj — > 7r~^ ir~ e~^ e~ {'y) decay using about 1.7 collected by the 
KLOE experiment at the DA$NE ^-factory. This corresponds to about 72 millions 
T] mesons produced in (p radiative decays. We have measured the branching ratio, 
inclusive of radiative effects, with 4% accuracy: BR{r] — > 7r^7r~e'''e~(7)) = (26.8 it 
0.9stat. i 0-7syst.) X 10~^. We have obtained the first measurement of the CP-odd 
vrvr — ee decay planes angular asymmetry, A(j, = (— 0.6± 2.5 stat. i 1.8 syst.) x 10~^. 

Key words: e~^e~ collisions, rare r] decays, CP violation 
PACS: 13.66.Bc, 14.40.Aq, ll.30.Er 



1 Introduction 

The decay of light pseudoscalar mesons, 7r°, r] and r]', proceeds via electromag- 
netic interaction and the radiative decays of t] and rj' to pions allow to probe 
their electromagnetic structure [1]. Conversion decays offer the possibility to 
measure precisely the virtual photon 4-momentum via the invariant mass of 
the e'^e~ pair. The branching ratio for internal conversion decay of the r) me- 
son, f] — s> 7r"*"7r~e^e~, has been computed with different approaches, but until 
recently both the theoretical and the experimental results were affected by 
large uncertainties. 

The first calculation, based on pure QED, is 40 years old: BR ~ 3 x 10^"^ [2]. 
The addition of vrvr interaction treated with the Vector Dominance Model gives 
BR = 3.6 X 10~^ with an error of about 10% [3], while an approach based 
on chiral perturbation theory that includes vector mesons gives (3.2 ± 0.3) x 
10~^ [4]. Recently, an approach based on the chiral effective Lagrangian in- 
cluding TTTT interactions has obtained a more precise result: BR = (2.99i'Iqq9) ^ 
10-4 [5]. 

The 7] — i> TT+TT-e^e- decay has been first observed by the CMD-2 experi- 
ment [6], giving BR = (3.7l^;8 ± 3.0) x lO'^, and has afterwards been con- 
firmed by the CELSIUS- WASA experiment [7,8]: BR = (4.3;?;^ ±0.4) x lO'^. 
The precision of these results does not allow to test different models. 
Recently, a possible CP violating mechanism (CPV), not directly related to 
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the most widely studied flavor changing neutral processes, has been pro- 
posed. This mechanism could induce interference between electric and mag- 
netic decay amplitudes. Such CPV effect could be tested in the decays of 
the pseudoscalar mesons by measuring the polarization of the virtual photon 
and would result in an asymmetry in the angle between the planes con- 
taining the e^e~ and the ii^Ti~ pairs in the meson rest frame, deflned as: 
= ilo^' %d<P-S:/2 %d<t^)/Uo" S# + /;/2 %d(t^)- This kind of asymmetry 
has been already studied [9-13] and observed [14,15] in the decay of the 
meson. In the rj decay this asymmetry is constrained, by experimental [16] 
and Standard Model [17] upper limits on the CP- violating decay r/ tt+tt^, 
to be, at most of O(10~^) and C(10~^^), respectively. However, as pointed 
out in [18,19], it is possible in case of sources of CPV beyond the Standard 
Model which do not contribute directly either to ex or to the neutron electric 
dipole moment, dn- In this case, is predicted to be up to C(10~^), a value 
reachable with the present statistics of ry mesons. 



2 The KLOE detector 

The KLOE experiment operates at DA$NE , the Frascati ^-factory. DA$NE 
is an e^e~ collider running at a center of mass energy of ~ 1020 MeV, the 
mass of the meson. Equal energy positron and electron beams collide at an 
angle of 7r-25 mrad, producing mesons nearly at rest. 

The detector consists of a large cylindrical Drift Chamber, surrounded by a 
lead-scintillating flber ElectroMagnetic Calorimeter (EMC). A superconduct- 
ing coil around the EMC provides a 0.52 T fleld. The drift chamber [20], 4 m 
in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense wires and 
37,746 aluminum fleld wires. The chamber shell is made of carbon flber-epoxy 
composite and the gas used is a 90% helium, 10% isobutane mixture. The 
position resolutions are a^y ~ 150 /im and cr^ ~ 2 mm. The momentum reso- 
lution is <j{p±)/p± ~ 0.4%. Vertices are reconstructed with a spatial resolution 
of ~ 3 mm. The calorimeter [21] is divided into a barrel and two endcaps, for 
a total of 88 modules, and covers 98% of the solid angle. The modules are 
read out at both ends by photo-multipliers, both in amplitude and time. The 
readout granularity is ~ (4.4 x 4.4) cm^, for a total of 2440 cells arranged 
in flve layers. The energy deposits are obtained from the signal amplitude 
while the arrival times and the particles positions are obtained from the time 
differences. Cells close in time and space are grouped into calorimeter clus- 
ters. The cluster energy E is the sum of the cell energies. The cluster time T 
and position R are energy-weighed averages. Energy and time resolutions are 
aE/E = 5.7%/^E (GeV) and at = 57 ps/y^E (GeV) © 100 ps, respectively. 
The trigger [22] uses both calorimeter and chamber information. In this anal- 
ysis the events are selected by the calorimeter trigger, requiring two energy 
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deposits with £" > 50 MeV for the barrel and E > 150 MeV for the end- 
caps. A cosmic veto rejects events with at least two energy deposits above 30 
MeV in the outermost calorimeter layer. Data are then analyzed by an event 
classification filter [23], which streams various categories of events in different 
output files. 



3 Event selection 



This analysis has been performed using 1,733 pb ^ from the 2004-2005 dataset, 
242 pb~^ from the 2006 off-peak (^i = 1000 MeV) data, 3,447 pb~^ of Monte 
Carlo (MC) simulating all (p decays, 50,506 pb^^ of signal Monte Carlo. Sig- 
nal MC has been generated according to the matrix element in [19], with 
BR = 4 X 10~^ and having Acf, = 0. All MC productions account for run by 
run variation of the main data-taking parameters such as background con- 
ditions, detector response and beams configuration. Data-MC corrections for 
calorimeter clusters and tracking efficiency, evaluated with radiative Bhabha 
events and ^ pvr samples respectively, have been applied. Effects of Final 
State Radiation (FSR) have been taken into account using the PHOTOS MC 
package [24,25]. This package simulates the emission of FSR photons by any of 
the decay products taking also into account possible interference between dif- 
ferent diagrams. We have inserted PHOTOS in our Monte Carlo at the event 
generation level, simulating FSR for the electrons, therefore our simulation 
fully accounts for radiative decays. 

At KLOE, rj mesons are produced together with a monochromatic recoil pho- 
ton through the radiative decay (p —>■ rj'-f [E^ = 363 MeV). In the considered 
data sample about 72 x 10^ ?7's are present. As first step of the analysis, a 
preselection is performed requiring at least four tracks (two positive and two 
negative) coming from the Interaction Point. The Fiducial Volume is defined 
by a cylinder having radius R = 4 cm and height h = 20 cm. For each charge, 
the two tracks with the highest momenta are selected. After track selection, 
one and only one neutral cluster, having energy E^i > 250 MeV and polar 
angle in the range (23° — 157°), is required. A cluster is defined neutral if it 
does not have any associated track and has a time compatible with the photon 
time of fiight. After preselection, the signal is about 1.4% of the sample. 
Mass assignment for each track is performed by either identifying a pion de- 
cay from a kink in the track, or using the Time Of Flight (TOF) of the 
particles. For each track associated to a calorimeter cluster, the quantity 
= ttrack — tduster in both elcctrou (Ate) and pion (Atjr) hypothesis is eval- 
uated; ttrack is defined as the length of the track divided by (5{m)c. The mass 
hypothesis minimizing At is then chosen. If two tracks of the same charge 
satisfy the same mass hypothesis, the minimum Atg identifies the electron. 
When only one of the same charge tracks is identified as pion (electron), the 
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other one is assumed to be electron (pion). For the remaining events the tracks 
with higher momentum are assumed to be pion. 

To improve the energy and momentum resolution, a kinematic fit is performed 
imposing the four-momentum conservation and the TOF of the photon. A very 
loose cut on the of the kinematic fit is applied in order to discard poorly 
reconstructed events. 



4 Background rejection 

Two sources of background have to be distinguished: 

(1) (f) background: 

this is mainly due to — ^ tt+tt^tt^ events (with vr^ Dalitz decay) and to 
— i> ?77 events either with rj tx^-k^-k^ (with 7r° Dalitz decay) or with 
7] 7r"^7r^7 (with photon conversion on the beam pipe). Note that this 
last background has the same signature as the signal. Background from 
— > KK is also present at the preselection level. 

(2) Continuum background: 

this is due to e^e~ —>■ e~^e~{'j) events with photon conversions, split tracks 
or interactions with some material in the region of DA$NE quadrupoles 
inside KLOE. Because of poor MC statistics, this background has been 
studied using off-peak data taken at y/s = 1 GeV, where decays are 
negligible. This sample has been properly normalized for luminosity and 
-y/s behaviour. 

A first background rejection is performed cutting on the sum of the absolute 
value of the momenta of the two particles having the highest momenta and 
opposite charge, s2p = \pipMax)\ + \piPMax)\y cutting on the sum of the 
absolute value of the momenta of the four selected tracks, s4p = X^i \Pi\- It is 
required that: (270 < s2p < 460) MeV and (450 < sAp < 600) MeV. 
To reject the events due to photon conversion on the beam pipe, we extrapolate 
backward the e~^e~ candidate tracks, down to the intersection with the BP, 
and compute there the invariant mass (Mge) and their distance [Dee)- A clear 
signal of photon conversion is visible in the Dee{BP)-Mee{BP) plane (figure 
1). We reject events having Mee{BP) < 15 MeV and D^eiBP) < 2.5 cm. The 
expected M^jree spectrum before and after this cut is shown in figure 2. The 
background events peaked at the r] invariant mass are significantly reduced 
by this cut. 

Finally, to remove continuum background from interactions with quadrupoles, 
we have defined the quantities {cos Of) and {cos6h) as the average polar angle 
of forward and backward particles identified as signal. Events having (cos 6'j) > 
0.85 and {cos 9b) < —0.85 are rejected. This cut affects neither the signal nor 
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Fig. 1. Dee vs Mee evaluated at the beam pipe for (p background (top panel) and MC 
signal (bottom panel). Events in the box Mee{BP) < 15 MeV Ci Dee{BP) < 2.5 cm 
are rejected. 
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Fig. 2. Spectrum of the T:~^'K~e'^e~ invariant mass after the cuts on the momenta 
(top panel) and after the cut to reject events with photon conversions (bottom panel) 
have been applied. The black histogram is the expected distribution, i.e. signal MC 
(dark grey), (j) background (light grey) and continuum background (white). The 
three samples have been normalized according to their luminosity. 



bacliground from decay events. 



5 Fit to the tt+tt e+e invariant mass spectrum and event count- 
ing 



In order to evaluate the background contribution, we perform a fit to the data 
distribution of the 7r"'"7r~e"'"e~ invariant mass after the cuts on the momenta. 
The fit is done on sidebands in order not to introduce correlations between 
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Fig. 3. vr+vr^e^e^ invariant mass spectrum on a wide range (left panel) and zoomed 
around the rj mass (right panel). The background scale factors have been obtained as 
described in section 5. Dots: data. The black histogram is the expected distribution, 
i.e. signal MC (dark grey), (j) background (light grey) and continuum background 
(white). 



signal and background. The ranges used are: [450, 520] Me^ U [570, 650] MeV. 
Upper and lower limits (450 and 650 MeV) have been chosen in order not 
to include in the fit tails from background distributions. The central range 
(520 and 570 MeV) is wide enough to well contain the tails of the signal 
distribution. Then the fit is performed using the background shapes only. 
The most precise background evaluation has been obtained fixing the off- 
peak data scale factor with luminosity at 7.14 ± 0.03 and fitting the MC 
background from (j) meson decay. The output of the fit is x^/dof = 32.5/30 
(P(x^) = 0.35), with a scale factor of 0.528±0.009, which is in good agreement 
with expectation from luminosity. The other possible approaches (fixing both 
the background scale factors with luminosity, leaving both free in the fit or 
fixing the decay and leaving free the continuum background shape) have 
been used for the evaluation of systematic uncertainties. The result of the fit 
is shown in figure 3 both in a wide M^^gg window and around the signal region. 
For the signal estimate we limit ourselves to the region [535, 555] MeV and 
perform the event counting after background subtraction: we find 1555 (368) 
signal (background) events, see table 1. Data-MC comparisons show a very 
good agreement for all considered variables. The most relevant distributions 
are reported in figure 4. 



6 Measurement of the BR(r7 — > tt+tt e'^e (7)) 

The branching ratio has been evaluated according to the formula: 

BRiv - 7r%-e+e-(7)) = ^^---'--^'^-(^^ 1 (1) 

-''r)7 >7r+7r~e+e~ (7) 
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Fig. 4. Spectra of opening angle (top) and invariant mass (bottom) for the e^e~ 
(left) and 7r"'"7r~ (right) pairs for events in the signal region. The background scale 
factors have been obtained as described in section 5. Dots: data. The black histogram 
is the expected distribution, i.e. signal MC (dark grey), (j) background (light grey) 
and continuum background (white). 

Table 1 

Event counting in the signal region. 



Contribution 


Counting 


Data 


1923.0 


Signal 


1555.0 


Total background 


368.0 


(j) background 


275.2 


Continuum background 


92.8 



where A''^_^7r+7r-e+e-(7) is the number of signal events and e,,^7r+^-e+e-(7) is 
the efficiency taken from MC. The number of —> 777 events, A'^^-y, has been 
obtained using the formula A^^^ = £ ■ cx^^^^, where £ is the integrated lumi- 
nosity and the cross section a^->n^ takes into account the meson line shape. 
a^^rn has been evaluated with 77 tt'^tt^tt^ events collected in the range 
1017 MeV < ^/s < 1022 MeV [26]. Inserting all the numbers summarized in 
table 2, we obtain the value: 

BR{t^ ^ 7r+7r-e+e-(7)) = (26.8 ± O.Qstat.) x 10"^ (2) 
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Table 2 

Summary of the numbers used in the master formula (1) for the branching ratio 
evaluation. 



BR inputs 


Values 


Number of events 


1555 ± 52 


Efficiency 


0.0803 ± 0.0004 


Luminosity 


(1733 ± 10) pb~^ 


e'^e~ ^ (/> — > r/7 cross section 


(41.7 ±0.6) nb 



where the error accounts for the event counting and the background subtrac- 
tion. 

The systematic uncertainties have been evaluated in the following way: 

• fixing with the luminosity or leaving free in the fit the background scale 
factors; 

• varying the sideband upper and lower limits in the ranges [400, 535] MeV 
and [555, 700] MeV respectively; 

• varying the binning of the histogram used for the fit from 1 MeV/bin to 5 
MeV/bin; 

• repeating the whole analysis chain after moving selection criteria on s2p, 
sAp, Dee{BP), Mee(5P), (cos^/), (cos ^fe) and M^^^e by ±1ct, ±2a's, iSa's 
around the reference value. The BR is then recomputed for all of these 
variations. The systematic uncertainty has been evaluated as the quadratic 
sum of RMS's obtained for each case. 

The uncertainty on A'^^^ has been also added to the systematics in the nor- 
malization term. The results of the systematics evaluation are summarized in 
table 3. The largest contributions are due to the normalization and to the cut 
on Mee(-B-P). Taking the total systematic error into account, the measurement 
of the branching ratio is: 

BR(^ ^ 7r+7r-e+e-(7)) = (26.8 ± 0.9stat. ± 0.7syst.) x lO'^ . (3) 



7 Decay plane asymmetry evaluation 

The decay plane asymmetry is calculated starting from the momenta of the 
four particles and is expressed as function of 0, the angle between the pion 
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Table 3 

Summary of the systematic uncertainties on the branching ratio. 



Source of uncertainty 


cj(BR) 


Pree/fixed scale factors 


0.18 X 10" 


-5 


Sidebands range 


0.05 X 10" 


-5 


Binning 


0.02 X 10" 


-5 


Analysis selection 


0.55 X 10" 


-5 


Normalization 


0.42 X 10" 


-5 


Total 


0.72 X 10" 


-5 




Fig. 5. Definition of the angle (p between the pion and electron decay planes, 
and the electron planes in the rj rest frame (figure 5): 

I -^sin (/)cos </)>0 -^sin (/) cos (/)<0 /^\ 

" Iv Tm ■ ^ ' 

' sin (/) cos (/i>0 ~r » sin (/)Cos 0<O 

The quantity sin cos is given by (nee x n.„T,)z{nee ■ nmr), where the n's are 
the unit normals to the electron and pion planes and z is the unit vector along 
the axis defined by the intersection of the two planes. The distribution of the 
sin cos variable in the signal region is shown in figure 6. We remind that 
the signal MC has been produced with = 0. 

While the analysis efficiency is completely fiat in the sin cos distribution, 
some distortion is introduced by the reconstruction, because of events with 
wrong mass assignment. The correction to this distortion has been evaluated 
by MC, fitting with a linear function the ratio between the generated and re- 
constructed sin cos distributions. The resulting slope is -0.032 ± 0.016. The 
use of higher polynomials does not improve the fit. The origin of this slope has 
been investigated by MC and it is completely due to the 14% of signal events 
with wrong particle identification. This because the particle identification al- 
gorithm forces the mass assignment in case of ambiguities without discarding 
events. The aim is to preserve the statistics, which completely dominates the 
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sini))cosit) 



Fig. 6. Distribution of the sin (p cos (p variable in the signal region. The background 
scale factors have been obtained as described in section 5. Dots: data. The black 
histogram is the expected distribution, i.e. signal MC (dark grey), cp background 
(light grey) and continuum background (white). 

asymmetry measurement. 

The asymmetry has been evaluated for the events in the 535 MeV < M^^ee < 
555 MeV mass region after background subtraction. After applying the cor- 
rection, we obtain: 

= (-0.6 ± 2.5 Stat. ± 1.8 syst.) X 10-2 (5) 
which is the first measurement of this asymmetry. 

As for the branching ratio, the systematic error has been evaluated repeating 
the whole analysis chain after varying selection criteria by ila, ±2o"'s and 
±3o"'s around the reference value and taking as uncertainty the quadratic sum 
of the resulting RMS's. The uncertainty due to the correction has been eval- 
uated varying its slope by ±la. The largest contribution is due to the cut on 
M^Tree while the contribution of the slope correction is 0.5 x 10"^. 
In order to check the distortion correction applied to A,/), we have defined 
a control sample having only events without ambiguities in particle identifi- 
cation. In this case the probability of wrong particle identification is almost 
zero and no distortions are observed in the MC sin0cos</) distribution. The 
fraction of this control sample in data and MC events is in good agreement 
(0.62 ± 0.02 and 0.64 ± 0.02 respectively), showing that our simulation repro- 
duces the real data well. The asymmetry evaluated with the control sample 
is in good agreement with our measurement but has a larger statistical error: 
= (-1.2 ±3.1 Stat.) X 10-2. 



8 Conclusions 

Using a sample of 1.7 pb^^ collected in the meson mass region, we have 
obtained a measurement of the rj 7T^7T~e~^e~{'-f) branching ratio with 4% 
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accuracy, ten times more precise than the previous best measurement [6-8]: 

BR{'q ^ 7r+7r-e+e-(7)) = (26.8 ± Q.^stat. ± ^-Tsyst) x IQ-^ . (6) 

Radiative events shghtly modify momentum distribution of the charged par- 
ticles and have been carefully considered in the efficiency evaluation. As a 
result, the measured branching ratio is fully radiation inclusive. 
Our measurement is about 2cr smaller than theoretical predictions [3-5], while 
it is in agreement (~ la) with the calculations of the ratio of the branching 
fractions BR{r] -k^ -k^ e~) / B Rij] 7r+7r~7) in references [2,5] when the 
recent CLEO measurement of BR{ri 7r"*"7r~7) [27] is used as normalization. 
The final sample of 1555 signal events allows us to perform the first measure- 
ment of the CP- violating asymmetry .4^, which is consistent with zero at the 
3% percent precision level: 

= (-0.6 ± 2.5 Stat. ± 1.8 syst.) X 10-2 . (7) 
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